Probing dark matter models with neutrinos from the Galactic center 



Arif Emrc Erkoca/ Mary Hall Reno,^ and Ina Sarccvic^''^ 

^Department of Physics, University of Arizona, Tucson, AZ 85721 
'^Department of Physics and Astronomy, University of Iowa, Iowa City, lA 5224-2 
^Department of Astronomy and Steward Observatory, University of Arizona, Tucson, AZ 85721 

We calculate the contained and upward muon and shower fluxes due to neutrinos produced via 
dark matter annihilation or decay in the Galactic center. We consider dark matter models in which 
the dark matter particle is a gravitino, a Kaluza-Klein particle and a particle in leptophilic models. 
The Navarro-Frenk- White profile for the dark matter density distribution in the Galaxy is used. 
We incorporate neutrino oscillations by assuming maximal mixing and parametrize our results for 
muon and shower distributions. The muon and shower event rates and the minimum observation 
times in order to reach 2cr detection significance are evaluated. We illustrate how observation times 
vary with the cone half angle chosen about the Galactic center, with the result that the optimum 
angles are about 10° and 50° for the muon events and shower events, respectively. We find that 
for the annihilating dark matter models such as the leptophilic and Kaluza-Klein models, upward 
and contained muon as well as showers are promising signals for dark matter detection in just a few 
years of observation, whereas for decaying dark matter models, the same observation times can only 
be reached with showers. We also illustrate for each model the parameter space probed with the 
2(7 signal detection in five years. We discuss how the shape of the parameter space probed change 
with significance and the observation time. 

PACS numbers: PACS; 95.35.+d, 14.60.Lm, 95.55.Vj, 95.85.Ry 



I. INTRODUCTION 

The identity of dark matter (DM), which constitutes 
23% of the total density of the Universe has been 
one of the most important open questions in astrophysics 
for more than seven decades Concerted theoretical 
and observational efforts are being made to detect the 
DM through photonic and Icptonic signals 0-113 ■ Exten- 
sions of the standard model of particle physics have been 
proposed to account for the observed anomalies in the 
PAMELA positron and FERMI electron plus positron 
data 0, H 11 i, [IS [iMl- These models also predict 
neutrino production via DM annihilation or decay (20l - 

Neutrino signals are complementary to photon or 
charged lepton signals from the Galactic center, or more 
generally from the Galactic halo. Neutrinos with energies 
of the order of the DM mass, < for < 10 TeV, 
propagate without being absorbed or deflected in transit 
toward the Earth. Neutrino telescopes have the capa- 
bilities to probe higher energies than the satellite based 
telescopes [2^, [l^. The approaches to DM searches, in- 
cluding the optimal angular coverage of the Galactic cen- 
ter or halo, are different depending on the details of the 
DM model and the location of the detector relative to 
the Galactic center. The strategies for uncovering neu- 
trino signals of DM decay or annihilation in the Galactic 
halo seen in underground or underwater detectors, are 
the subject of this paper. 

We calculate the neutrino induced contained and up- 
ward muon flux, hadronic shower flux and the muon and 
shower event rates for different DM models which are 
considered to explain data excesses in other indirect DM 
searches. In this paper, we study the annihilation of the 



lightest Kaluza-Klein particle [13, [s^l, leptophilic DM 
particle annihilation [3, H^, two-body and three-body 
decay channels of gravitino, as well as the decay of the 
leptophihc DM [ia, [TMl m, Hi- We show that for 
each model, the shape of the muon and shower fluxes 
differs significantly from the shape of the neutrino fiuxes 
at production due to the smearing produced by neutrino 
interactions and muon propagation as discussed in Rcf. 
[31I [3^ . We also calculate the detection significances of 
the DM signals at IceCube/DeepCore [l^ detector and 
compare the different DM models in terms of the energy 
spectra and the total counts of the muon and shower 
events. We show the parameter space of different DM 
models that is being probed when we require 2a muon 
or shower signal within five years of observation. 



II. MODELS FOR DARK MATTER 

The ingredients for theoretical predictions of particle 
fluxes from DM annihilation or decay include a model 
for the DM distribution, a particle physics model for the 
DM particle couplings to standard model particles, and 
standard model physics processes for the resulting pro- 
duced particles. For the DM distribution in the galaxy, 
we use the Navarro-Frenk- White (NFW)[3^ profile as a 
typical realistic DM density profile. The expressions for 
the neutrino fluxes and its dependence on the DM profile 
are presented in Appendix A. In case of DM decay (an- 
nihilation), the neutrino spectrum has linear (quadratic) 
dependence on the DM density. 

The particle physics models on which we focus in this 
study consist of the leptophilic, Kaluza-Klein and grav- 
itino dark matters. Either thermal averaged annihilation 
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cross section times velocity {av) or a decay time r spe- 
cific to the model is required. Characteristically for an- 
nihilation, the required {av) is larger [l^, [l^ [22, IH, [l^l 
than the value required for a thermal relic abundance 
[11]: (ctw)o = 3 X 10^^^ cm^s~^. Following the current 
convention, we write 



{av) = B {av)o 



(1) 



with a boost factor B. There arc theoretical evaluations 
of the boost factor [sH, however, we treat the boost fac- 
tor as a phenomenological parameter in this paper. To 
explain the lepton excesses, some models have constraints 
on the boost factor as a function of DM mass • 

In Icptophilic DM models [3, HB] explaining the 
PAMELA positron excess, the DM annihilation or decay 
must proceed dominantly to leptons in order to avoid 
the overproduction of antiprotons. Moreover, according 
to the FERMI data, the direct production of electrons 
must be suppressed with respect to the production of 
electrons (and positrons) as secondaries. It was shown 
[IB] that the leptophilic DM with mass (m^) in the range 
between 150 GeV and a few TeV, which annihilates or 
decays into r's or /i's can fit the PAMELA Q and Fermi 
[^ data as well as the HESS high energy photon data [1] . 
The best fit parameters for the boost factor [B) and the 
decay time (r) which determine the overall normaliza- 
tions, for the specific case involving muons from annihi- 



lation [xx 
are given by [l6| 

B 



) or decay (x 



), respectively. 



43l77ix - 38.9 
2.29 X 10^6 sec 



sec 



^ Br X 10^ 



(2) 



for in TeV. The annihilation channel into tau pairs 
is less favored by the data [l6| . 

Some Kaluza-Klein models can provide a DM candi- 
date which gives the correct relic density [s^. To ac- 
count for the HESS results [^, the lightest Kaluza-Klein 
particle (LKP) would have a mass of the order of a TeV 
[14| . The LKP is also assumed to be neutral and non- 
baryonic. In this model, the particle couplings are fixed 
such that LKP pairs annihilate into quark pairs (35%), 
charged lepton pairs (59%), neutrinos (4%), gauge bosons 
(1.5%) and higgs bosons (0.5%) [11 [s^]. 

The first DM candidate proposed in the context of su- 
persymmetry is the gravitino ('03/2) which would be the 
lightest supersymmetric particle (LSP). The gravitino is 
the superpartner of the graviton. With the existence of 
small R-parity breaking to allow the LSP to decay, the 
gravitino decays into standard model particles. The de- 
cay rate of the gravitino in this scenario is so small that 
it can have a sufficiently long lifetime for the correct DM 
relic density today. 

In order to account for the observed anomalous 
positron excess in the PAMELA data and positron plus 



electron excess in the FERMI data, the lifetime of the 
gravitino DM is constrained to be of the order of ~ 10^® 
seconds and its mass to be in the range between few 100 
GeV and few TeV 17[. To explain the data, the three- 
body gravitino decay mode (V'3/2 ~^ was consid- 
ered [13] . We use the parameters of this model to explore 
neutrino signals from gravitino decay. For illustration, in 
addition to three-body decay, we also consider the two- 
body gravitino decay modes (7/'3/2 — > {W^l"^ , Zi',-fv)) 
assuming the same lifetime and mass as for the three- 
body decay, and with the branching fractions given in 
Table I. 
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TABLE I: Branching fractions for the two-body gravitino 
decay into different R-parity violating channels for different 
masses UM- 



Particle/mode 



mass Bt or B 



'3/2 



400 GeV Br=2.3 
ips,/2^ {Wl,Zu,'~iv) 400 GeV Br=2.3 
X ^^^^J.~ 2 TeV Br=2.9 
5(1)5(1' ^ {qq,l+l-,W+W-,ZZ,uu) 800 GeV B = 200 
XX 1 TeV B = 400 

TABLE II: Model parameters characterizing fits to explain 
FERMI and PAMELA anomalies used as examples in this 
paper. 

Selected DM model parameters are shown in Table HIl 
For each of the DM models considered, the decay dis- 
tribution of the produced particles to neutrinos in case 
of DM annihilation, or the gravitino decay distribution 
to neutrinos, enters into the calculation of the neutrino 
fluxes that arrive at Earth. For annihilation directly to 
neutrinos, the energy distribution of each neutrino is a 
delta function in energy, with the energy equal to the DM 
mass. This case has been well studied in the literature 
[21I [sil . [321 . [H]. Here, we look at the secondary neu- 
trinos. Fig. [1] shows neutrino spectra, plotted in terms 
oi X = / Ejj^rnax whcrc E^^max = for annihilating 
DM and E,y^max ~ ^x/"^ decaying DM models. The 
curves in the figure are normalized to count the number 
of neutrinos, and in the case of the Ziy final state, the 
fraction of Z decays to neutrinos. The muon neutrino 
spectra in the figure should be multiplied by the branch- 
ing fraction for a specific decay channel in a given model. 
Analytic expressions for the neutrino spectra are given 
in Appendix B. 
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FIG. 1: Muon neutrino (v^) spectra in terms of x = 
Ev / Ejj^rnax from the three-body decay of gravitino (dot-dash- 
dashed line), from the decay of r (dashed line), Z boson (solid 
line) and from one of the two-body decay channels of gravitino 
(i/) Zv) for which Breit-Wigner distribution is used. The 
distributions should be multiplied by the branching fractions, 
and oscillations should be taken into account for the flux of 
neutrinos at Earth. 



FIG. 2: Muon neutrino [u^j) fluxes from the annihilation of the 
Kaluza-Klein (dotted line), leptophilic (dashed line), and the 
decay of leptophilic (dash-dot-dashed line), three-body decay 
(dot-dashed line) and two-body decay (dot-dot-dashed line) of 
gravitino DM particles. Neutrino oscillations have been taken 
into account. The angle-averaged atmospheric muon neutrino 
flux at the surface of the Earth is also presented (solid line). 
The corresponding values of the parameters for each model 
are shown as well. 



III. NEUTRINO FLUX 

The neutrino flux at Earth can be evaluated using the 
neutrino flux expressions given by Eqs. (|A1[) and (|A2[) 
in Appendix A and the neutrino spectra given in Ap- 
pendix B with taking the neutrino oscillation effects into 
account. In Fig. [21 we show the muon neutrino flux 
at Earth for three decay channels and two annihilation 
channels. 

In some DM models all three flavors of neutrinos can 
be generated by DM annihilation or decay, implying the 
flavor ratio at the production to be : : Vt is 1:1:1. 
This ratio remains unchanged with oscillation. This is 
the case for the gravitino decay and Kaluza-Klein DM 
annihilation. However, in case of the leptophilic DM 
model, in which x ~^ the initial neutrino flavor 

ratio is 1 : 1 : which becomes 1 : 0.5 : 0.5 as neutrinos 
travel astrophysical distances. We take this oscillation ef- 
fect into account when we evaluate muon neutrino fluxes 
presented in Fig. [5] and when we evaluate muon event 
rates below. 

Fig. [2] shows that with the exception of the gravitino 
decays, the distributions of neutrinos have very weak en- 
ergy dependence. The two-body gravitino decay gives a 
spiked feature at the kinematic limit in neutrino energy. 
The relative normalizations of the DM curves comes from 
different DM lifetimes or boost factors. 

Also shown in Fig. [2]is the angle-averaged atmospheric 



muon neutrino flux at the surface of the Earth. It is char- 
acterized by an approximate formula [sil [37| , (in units 
of GeV~^km~^yr~^sr~^) 



= NoE,-^-^\ -r_ln(l + 6S,) + 



ln(l + eE^) 



(3) 



where the parameters in the formula are listed in Ta- 
ble mil These same parameters appear in the angle- 
dependent atmospheric neutrino flux for zenith angle 9, 



dEydO. 



NoE,--^-' 



1 + bEyCosO 1 -|- eE^cost 



(4) 



This formula does not account for the prompt neutrino 
flux [ssj . however, for the energy range of interest, the 
prompt atmospheric neutrino flux is negligible. 

The angle-averaged atmospheric neutrino flux is a good 
approximation. In Fig. [31 we show the angle-averaged 
flux from Eq. (4) and the flux from Eq. (5) with 6 = 60° 
and the integrated flux measured by the AMANDA-II 
detector from Ref . [s^ . The angle-averaged flux is a bit 
larger than the flux at 60°, so at least for 6 less than 60°, 
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1.74 


a 


0.018 
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0.024 GeV"^ 
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0.0069 
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0.00139 GeV"' 


No 


1.95 X lO^'' for u 
1.35 X 10^'^ for 77. 



TABLE III: Parameters for the atmospheric i'^ and ly^ fluxes 
^iven by Eqs. Q and Q, in units of GeV~^km~^yr~^sr~^ 

using the angle-averaged atmospheric flux gives a small 
overestimate of the atmospheric background. 



ing the boost factor. A comparison of the upward-going 
muon rate, where the target volume is enhanced by the 
muon range at high energies, and the contained rate of 
muon production by neutrinos in the detector, is a useful 
exercise. 

For the IceCube/DeepCore detector, the Galactic cen- 
ter is above the horizon, so the upward muon rate of DM 
produced neutrinos is from the Galactic halo in a direc- 
tion pointing away from the Galactic center. We consider 
this possibility as well. 

Showers, either electromagnetic or hadronic, are pro- 
duced by neutrinos. We look at the optimization for these 
as well as a function of cone half angle, but we note that 
the current capabilities for shower angular resolution are 
somewhat limited, on the order of 50° [40l. |4]|. 



A. Muons 



t AMANDA-II 
--■ ATM 8 = 60° 
— ATM angle-averaged 




FIG. 3: Angle-averaged atmospheric muon neutrino (v^+V^) 
flux (solid line) and the atmospheric flux for fixed 9 = 60° 
(dashed line) compared with the angle-averaged (fu+Vu) &ux 
from AMANDA-II muon neutrino flux measurements [3911. 



IV. NEUTRINO SIGNALS 

For each DM candidate and model, there are several 
signals to pursue in underground detectors. One possibil- 
ity is to measure or constrain the rate of muons produced 
by muon neutrinos, over and above the expected atmo- 
spheric background rate. High energy muons point es- 
sentially in the same direction as the incident neutrino, 
and the angular resolution of high energy muon tracks 
is quite good. With good enough energy and angular 
resolution, and a large enough target volume, one looks 
for neutrinos coming directly from DM annihilation in 
the Galactic center, however, the target volume may be 
a limitation for constraining model parameters includ- 



The neutrinos coming from the Galactic center and 
Galactic halo can produce muons through charged cur- 
rent interactions in the detector (contained muons) . The 
flux is given by 



dE, 



dcTZ{E,,E^ 



d<i)u \ Nap 
dE,, 



dE„, 



+ {p^7l) (5) 



where Na = 6.022 x 10^^ is Avogadro's number, p is 
the density of the medium, Emax = "nix ^'^^ annihilation 



and Err 



^/2 for decay. The differential cross sec- 



tions da^'^ /dE^ are the weak scattering charged-current 
cross sections for neutrino and antineutrino scattering 
with protons and neutrons 
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We evaluate the muon flux from neutrino charged- 
current interactions in the detector, when neutrinos are 
produced in DM annihilation or DM decay. In Fig. H] 
we show muon fluxes for the case when the DM parti- 
cle is a gravitino, a Kaluza-Klein particle and for a lep- 
tophilic model in which DM annihilation or decay pro- 
duces fj,'^ p~ , for the model parameters listed in Table II. 
We take the cone half angle around the Galactic center 

to be ^max = 1°. 

In case of the gravitino DM decay and for Kaluza-Klein 
DM annihilation, there are discontinuities in the slopes 
at the highest muon energies coming from the superpo- 
sition of the direct neutrino production (dot-dot-dashed 
line for DM decay and dotted line for the Kaluza-Klein 
annihilation). The direct neutrino production, xx ~^ 
is the "golden channel" for DM detection because in this 
case the muon flux is increasing with energy, and it peaks 
at Ey, ~ [3^ . 

As noted in Sec. II, the parameters used for DM 
masses, boost factors and lifetimes are characteristic of 
those that were shown to describe PAMELA, Fermi/LAT 
and HESS data [H, [H, [13 ■ Changing the value of the 
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FIG. 4: Muon flux for the contained events for gravitino decay 
(dot-dashed and dot-dot-dashed lines), Kaluza-Klein annihi- 
lation (dotted line), leptophilic model (dashed line for anni- 
hilation and dash-dash-dotted line for decay) compared with 
the atmospheric background (solid line), for the case when 
^max = 1°- Model parameters are given in Table II. 



boost factor or the lifetime affects only the overall nor- 
malization of the muon flux. We find that for this choice 
of the parameters, DM signals in leptophilic model ex- 
ceed the atmospheric background for E"^ > 175 GeV, 
while for the Kaluza-Klein DM model the signal is above 
the background for E"^ > 275 GeV. In both cases, the 
signal cuts off when = m^. 

We consider the effect on the muon flux shape when 
we change the parameters, for example for the leptophilic 
model. In Fig. [S]we show the contained muon flux from 
DM annihilation and decay in a leptophilic model for 
different values of the parameters B, r and m^^ which are 
constrained to satisfy Eq. ^ to describe the data fl6| . 
The decays (lower thin lines) have lower fluxes than the 
annihilations (upper thick lines), even though the shapes 
are similar. For leptophilic models, one cannot enhance 
the signal rate by increasing B or decreasing t with m-f^ 
fixed if the Fermi and PAMELA data are explained by 
the model. 

Contained muons, produced by neutrino interactions 
in the detector, make up one set of muon signals. Muons 
can also be produced in neutrino interaction in the rock 
below the detector. Muons produced with energy i?^, 
interact with the medium and finally reach the detector 
with energy E^. The effective volume of the detector is 
enhanced by the muon range at high energies. We denote 
these events as upward muon events. 

The muon range in the rock, i?p(i?^, E^), depends on 
the initial muon energy i?^, the final energy and 
the parameters a and /S which characterize muon en- 



FIG. 5: Contained muon flux from the annihilation, XX ~^ 
fJ-^ l-L~ (thick lines) and the decay, x ~^ (thin lines) 

processes. The relation between the boost factor and m^, 
and between the lifetime and m-^ are given by Eq. ([2]). We 

take 6max = 1°. 



ergy loss. Numerically, a ~ 2 x 10^'^ GeVcm^/g ac- 
counts for the ionization energy loss and /3 ~ 3.0 x 10"^ 
cm^/g for the brcmsstrahlung, pair production and pho- 
tonuclear interactions. The range is then approximated 
by R^iEl^E^) = \n[iEl + a/(3)/iE^ + a//3)]//3p. For 
muon transit through the rock, the muon range is 1 km 
for Ei ~ 1 TeV. 

Taking into account the energy losses, the final muon 
flux at the position of the detector can be written as [3l[ , 



dE„ 



dE^ 



dE„ 



2 V dE^ 



(6) 



where Psurv{E^^, E ^) is the survival probability for a 
muon with initial energy E^ to reach final energy E^. 
A detailed derivation of Eq. (6) can be found in Rcf. 
[sH ]. To first approximation, Psurv — 1 since at high 
energies, the muon has a long decay length. The muon 
energy at the production point is related to the muon 
energy a distance z from that point by 



Eliz) c e^P^E, + {e^"^ ~ 1) 



/3 



(7) 



We show in Fig. [51 the upward muon fiux for a generic 
northern hemisphere detector, looking down through the 
Earth with a cone half angle of ^^Inax = 1° around the 
Galactic center. The muon fluxes are smoothed rela- 
tive to Fig. m as a consequence of the energy loss. For 
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FIG. 6: Upward muon flux for the anniliilating and decaying 
DM models from Table II. We take e^ax = 1°. 

DM particles with masses of order 1 TeV, the upward 
muon flux is larger at low energies than for the contained 
muons because the muon range is larger than 1 km, effec- 
tively enhancing the volume of the kilometer-size detec- 
tor. When = 400 GeV, the energies of the produced 
muons are such that the muon range is less than 1 km, 
which is the size of the detector for which the contained 
muon flux was calculated. If the depth of the detector is 
500 m, the contained and upward muon fluxes for the DM 
mass of 400 GeV would be approximately equal at low 
energies, although the contained muon flux would have 
a little harder spectrum. This is direct consequence of 
muon range dependence on the DM mass. For example, 
muon with initial energy of 400 GeV (1 TeV) has a range 
of 500m (1 km). 

The contained muon event rate, iV^*(m^), is obtained 
by integrating muon flux folded with the effective volume 
of the detector, Vcs, i.e. 

Ki^x) = r^'" ^VME,)dE, (8) 

where dcj)')^ / dE^ is given in Eq. ((Sj and E^^ is the muon 
detector threshold, typically 10-100 GeV for deep ice 
or water detectors [28|. In our calculations, we choose 
E^J^ = 50 GcV. We also consider an energy independent 
IceCube/DecpCore effective volume, Vas — 0.04 km'^, for 
the contained muon events [l^, . 

Similarly, the upward muon event rate, is obtained by 

K^i^x) = / ^A,s{E^)dE, (9) 

where d(j)^P/dE^ is given by Eq. A^e is the angle- 
averaged muon effective area for which we assume A^s = 
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1 km2. 

The event rates for contained and upward muons for 
a cone half angle of ^^max = 1° are shown in Table IIVI 
for the DM models shown in Figs. H] and [Hand listed in 
Table II. We also obtain the number of years required for 
the rate of signal events s and background events b to 
satisfy the condition 

Vs + b 

From Table IIVI we note that for the parameters that we 
considered, only the Kaluza-Klein DM and leptophilic 
annihilation models have a reasonable chance of detection 
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TABLE IV: Event rates per year and the time required to 
reach 2a detection significance for the upward and the con- 
tained muons (^) for a cone half angle of Smax = 1°. Results 
for different DM models are obtained by taking AbH = Ikm^ 
and VcB = 0.04 km"* for the upward and contained muon 
events, respectively. 

The model parameters such as DM masses, annihila- 
tion cross sections and decay times that we consider are 
introduced to explain some indirect DM searches as ex- 
plained in the previous section. However, it is also pos- 
sible that the signals that have been observed (^-Q in 
these searches have no DM origin. Then, and B or 
T can be varied independently. In terms of the neutrino 
signals, the dependence of the signals on the annihilation 
cross sections or on the decay times is trivial since these 
parameters affect only the overall normalization. The 
dependence on DM mass is not that straightforward. 

In order to see the dependence of the signals on DM 
mass, we set the values of the boost factor and the decay 
times to those in Table II and calculate the event rates 
as a function of DM mass for each model. We present 
our results for the contained muon event rates in Fig. [7] 
and for the upward muons in Fig. [S] The solid line in 
each figure corresponds to the muon background due to 
the atmospheric neutrinos. 

From Fig. [7] we note that the contained muon events 
rates for annihilating DM models decrease with m^. On 
the other hand the event rates for the decaying DM mod- 
els increase slowly with m^, for < 1 TeV and for 
> 1 TeV, they become almost independent of m^. 
The m-^ dependence of the contained muon event rates is 
mainly due to the iti~'^ (i^^x^) dependence in the neu- 
trino flux for DM annihilation (DM decay) combined 
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FIG. 7: Contained muon rates as a function of DM mass for 
tlie models presented in Table II. We take the muon detector 
threshold to be Ejl' = 50 GeV and the cone half angle to be 

n 1 o 

"max — ^ • 
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FIG. 8: Upward muon rates as a function of DM mass for the 
DM models presented in Table II. We take the muon detector 
threshold to be Ejl^ = 50 GeV and the cone half angle to be 

dmax — 1 • 



with the upper limit of integration dependence on m^. 
For DM masses in the range Ejj^ < < 400 GeV, 
where E*^ = 50 GeV, the integration region is sensitive 

to the value of the DM mass, while for to^ ^ Ej^ there 
is only weak dependence on m^. These combined effects 
are responsible for the observed m-y dependence of the 



The dependence of the upward muon rates is shown 
in Fig. [HI We find that the event rates for decaying DM 
models increase with while for annihilating DM mod- 
els there is almost no dependence for a wide range of 
DM masses. In contrast to the contained muon rates, for 
upward muons there is additional that is present in 
the muon range. As we increase the value of DM mass, 
the effective volume which depends on the muon range 
in rock becomes larger. 

The upward muon rates for a decaying DM particle 
have steeper increase with increasing DM mass than for 
contained muon rates, because of the energy dependent 



effective volume which increases with 



when com- 



pared to the case for the contained muon events. 

In Fig. [HI we present results for DM annihilation cross 
section required for a given DM mass in order to reach 2a 
detection significance in five years of observation within 
^'max = 1° for Kaluza-Klein (solid lines) and annihilating 
leptophilic (dashed lines) models. From Fig. [Tjwe note 
that the contained muon event rates decrease with for 
a fixed annihilation cross section (i.e. fixed boost factor). 
Therefore, in order to have the same detection signifi- 
cance for each DM mass, DM annihilation cross section 
needs to increase with m-^ as shown in Fig. [5] However, 
for the upward muons, the event rates increase with 
for < 1 TeV and exhibit a slight decrease for higher 
DM masses for a fixed annihilation cross section. Thus, 
in order to have the same significance independent of the 
DM mass for the upward muon events the DM annihila- 
tion cross section has to decrease with for < 1 
TeV and increase for > 1 TeV as seen in Fig. [9l If 
there is no signal detected at 2a level in five years, the 
parameter space above each curve is excluded at that sig- 
nificance level. Our results also indicate that the upward 
muons are more promising than the contained muons in 
constraining the model parameters. Increasing the ob- 
servation time would result in larger excluded parameter 
space. 

Similar to the models for the annihilating DM, we eval- 
uate the parameter space for the decaying DM models. 
In Fig. [TUj we show the decay time as a function of the 
DM mass for 0max = 1° that is needed in order to reach 
2a detection significance with five year observation pe- 
riod. For a wide range of DM masses, the contained 
muon event rates have weak dependence on m^, while 
the upward muon event rates show a steep increase with 
increasing to^ as can be seen from Figs. [7] and [8] This 
implies that we need longer decay time for the upward 
muon events than for the contained muon events to reach 
the same detection significance for a five year observation 
time while the decay time has almost no dependence on 
the DM mass for the contained muon events. The param- 
eter space below each curve corresponds to the exclusion 
region at 2a level after five years of no signal detection. 
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FIG. 9: Annihilation cross section versus DM mass for upward 
and for contained muon events to reach 2cr detection signif- 
icance with five years of observation for the case of Kaluza- 
Klein (solid lines) and leptophilic (dashed lines) models. We 
take Smax = 1° and = 50 GeV. 



B. Optimal angles and muon signals 

As noted earlier, size of the cone (cone half angle value) 
that is too small can restrict the number of events too 
much. In Fig. 111! we present contour plots similar to 
those in Fig. [10] but for ^^ax = 10°. With increasing 
Smax from 1° to 10°, the J factor ((Ji)oAi7) entering 
into the neutrino flux increases, resulting in stronger re- 
striction on the DM lifetime. 

For more general use in discussing the dependence on 
cone half angles, the values of (Ji)oAi7 and (J2)oAf2 
given in Table VI in Appendix A can be combined with 
the parametric forms for the muon fluxes given in Ap- 
pendix C to calculate event rates for different cone half 
angles around the Galactic center. From Table VI we 
note that the cone size dependence of the J factors is 
quite different for the case of DM annihilation than for 
decaying DM. 

The difference between DM density contributions to 
DM annihilation and decay can also be seen in Fig. [20] in 
Appendix A, where we present ( J„)nAJ7 factors for DM 
annihilation {n = 2) and DM decay {n = 1) evaluated 
for a cone wedge between ^max — 1° and ^max around the 
Galactic center. We note that moving the wedge away 
from the galactic center leads to a significant reduction of 
the signal from annihilating DM due to the dependence 
of (J2)nAri on square of the density (see Eq. (|A3p ). We 
find that the signal is reduced by a factor of 17 when the 
wedge is moved from 1° to 90° off the galactic center. On 
the other hand, in a wedge between 50° and 90° around 
the Galactic center, the value of (J2)f2Ari is about 5.6, 



FIG. 10: Decay time versus DM mass for upward and for 
contained muon events to reach 2a detection significance with 
five years of observation for the case of decaying DM models: 
gravitino three-body (solid lines), gravitino two-body (dashed 
lines) and leptophilic (dotted lines). We take Smax = 1° and 
Ej!" = 50 GeV. 



a value close to what one obtains for a cone centered at 
the Galactic center with Smax = 5°. 

In contrast, for a decaying DM, the signal from a wedge 
between ^max— 1° and ^max around the Galactic center in- 
creases with 6'inax for 6'max < 30° and decreases slowly for 
higher 6'inax- For example. {Ji{l°))aAfl = 0.018 whereas 
(Ji(30°) - Ji(29°))oA17 0.2 which is only a factor of 
2 higher than (Ji(90°) - Ji(89°))nAf} = 0.1. This is a 
consequence of the dependence of the signal for a decay- 
ing DM particle on density which is one power less when 
compared to an annihilating DM particle. An increase 
in the size of the wedge (i.e., the volume of the source 
region) can result in an enhancement in the signal even 
if the chosen wedge is away from the Galactic center. 

In what follows we consider cone half angles larger than 
6'max = 1°. In Fig. [T2] the detection times required for 
the signal s to satisfy Eq. (fTU]) for the Kaluza-Klein and 
leptophilic annihilation cases are shown as a function of 
cone half angle. The optimal angle depends on the model 
and is of the order of a few degrees. The observation time 
is almost independent on the cone half angle in the range 
of a few degrees to about 10 degrees, but it increases for 
larger cone half angles where the signal increases slower 
than the background. 

As shown in Fig. [T2]for the decaying DM models, the 
observation time to reach the same significance decreases 
with the cone half angle ^max, for ^max < 50°, and it 
increases only slowly for the higher cone sizes, as the 
signal in directions away from the Galactic center are 
relatively more important. 
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FIG. 11: Decay time versus DM mass for both upward and 
contained muon events to reach 2a detection significance after 
five years of observation for the case of decaying DM models: 
gravitino three-body (solid lines), gravitino two-body (dashed 
lines) and leptophilic (dotted lines). We take Smax = 10° and 
E'J' = 50 GeV. 



FIG. 12: The observation time t versus cone half angle Smax 
centered at the Galactic center to reach 2a detection signif- 
icance for the upward and contained muon events produced 
by the muon neutrinos originated from the annihilation of 
DM particles. Kaluza-Klein (solid lines) and leptophilic where 
XX A'^M" (dashed lines) models are considered. We take 
Ef = 50 GeV. 



We find that it takes, in general, shorter amount of 
time for the upward muon events than the contained 
muon ones to reach the desired detection significance. It 
is mostly because of the small cfFcctive volume of the de- 
tector for contained events which results in lower event 
rates. The background is also small in this situation, 
however, still longer observation times are required due to 
low statistics. Among all the decaying DM models that 
we consider, the leptophilic model (x —> /i+/i^) seems 
most promising for detection of DM signal at 2a level for 
^max > 20° via upward muon events within a few years 
of observation. 

C. Hadronic and Electromagnetic Showers 

In addition to muons, the showers (hadronic and/or 
electromagnetic) produced in neutrino interactions could 
be used as signals of DM. The shower flux is given by 

dEsh Je^, \dEj 2 

(dal{E,,E,~E,h) , , A 

+ i^-^^). (11) 

where da'''" /dEgh is the differential cross section for 
showers produced in neutrino and antincutrino charged- 



current and neutral-current interactions with protons 
and neutrons. 

In Fig. 1141 we show shower flux for the case when 
the DM particle is a Kaluza-Klein particle, a gravitino 
and for a leptophilic model, for the model parameters 
listed in Table II. We take the cone half angle around 
the Galactic center to be 6'max = 50°. The shapes of the 
shower fluxes are similar to the contained muon fluxes 
presented in Fig. 4. The "golden channel" , XX ~^ 
in the case of contained muons, contributes to the shower 
events at the highest kinematically allowed energy for 
the case when the DM particle is a Kaluza-Klein particle 
(dotted line) and when it is a gravitino (dot-dot-dashed 
line). While for the case of contained muons (Fig. 4), 
DM signal in models with annihilating DM (Kaluza-Klein 
and leptophilic) exceeded the atmospheric background 
at high energies, shower signals for the same models are 
below the background. In addition, shower signals for the 
models in which DM annihilate are comparable to those 
in which DM decays. This is mainly due to our choice of 
large 6'max for the showers (6'niax = 50°). This can be seen 
by comparing (J2)si (for DM annihilation) and ( Ji)si (for 
DM decay) given in Appendix A (Eq. (|A3P ) for different 
values of the cone size. For example, {J2)nl {Ji)n — 75 
for 6'max 1° whereas (J2)n/(Ji)o = 3 for 6'max = 50°. 
Furthermore, shower fluxes extend to higher energies for 
the models in which DM annihilates than the case of DM 
decay, due to the kinematic constraints. 

We evaluated the shower rates by integrating Eq. PT|) 
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FIG. 13: The observation time t versus cone half angle 0max 
centered at the Galactic center to reach 2a detection signif- 
icance for the upward and contained muon events produced 
by the muon neutrinos originated from the decaying DM par- 
ticles. Decay of leptophilic where x ~^ M^A*" (dashed lines), 
two-body decay (solid lines) and three-body decay (dotted 
lines) of gravitino are considered. We take Ej^ — 50 GeV. 



FIG. 14: Flux for the shower events for two-body (dot-dot- 
dashed line) and three-body (dot-dashed line) gravitino decay, 
decaying leptophilic model (dash-dash-dotted line), Kaluza- 
Klein annihilation (dotted line) and annihilating leptophilic 
model (dashed line) compared with the atmospheric back- 
ground (solid line), for the case when Smax ~ 50°. Model 
parameters are given in Table II. 



over the shovirer energies 

N,him^) = [ ""^ ^V,s{E,h)dEsh, (12) 
J^th atjsh 

sh 

where T4ff (-Es/i) is the effective volume of the detector 
for measuring showers. 

In Table |Vl we present the shower rates for a cone size 
of 6',„ax = 50° for the DM models shown in Fig. [H For 
the detector volume, we use an energy independent Ice- 
Cube/DeepCore volume whose value is T4ff = 0.02 km^ 
for the showers [H H^, 113. The detector threshold is 
again set to be = 50 GeV. Although the number 
of events for all models is relatively small compared to 
the atmospheric background, the number of years needed 
for detection of the 2a effect is encouraging for the lep- 
tophilic model as seen in Table [V] We note that showers 
seem to be the best way to look for the signal of decaying 
DM. 

The shower event rates depend on the DM decay time 
and the boost factor as an overall normalization while the 
dependence on DM mass is non-trivial. In order to illus- 
trate shower event rate dependence on the DM mass, we 
take fixed boost factor or fixed decay time used in Fig. 

for the annihilating (decaying) DM models and we 
obtain shower rates for different m^. From Fig. 1171 we 
note that the shower event rate dependence on has 
similar behavior to to the case of contained muon events 
(see Fig [7]) , due to the similar flux of showers and muons 
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TABLE V: Shower event rates per year and the time required 
to reach 2a detection significance for the hadronic showers. 
We consider DM signals from different models described in 
the text. We take the effective volume to be VeS — 0.02km^. 



in neutrino charged-current and neutral-current interac- 
tions. For the case when the DM is a Kaluza-Klein parti- 
cle or in the leptophilic (xx ~^ model the shower 
rates decrease with the DM mass while for the decaying 
DM models it is slowly increasing with for < 1 
TeV and become almost independent for higher val- 
ues of m^. For the shower events, the differential weak 
scattering cross section in Eq. pT|) becomes higher for 
low shower energies so that the rates for the showers tend 
to increase more than those for the muons as the energy 
decreases. Consequently, the dependence on the choice of 
the detector threshold becomes more significant for the 
DM masses close to the threshold energy. We also note 
that within the cone half angle of 50° and with the chosen 



11 



10 F 




e (degrees) 

max ' 



FIG. 15: The time t versus cone half angle 0max centered 
at the Galactic center to reach 2a detection significance for 
the hadronic showers produced by the neutrinos from XX ~^ 
fi'^ IJ.~ (dashed line) and from the Kaluza-Klein DM particle 
(solid line). We take Ef^ = 50 GeV. 



model parameters the rates for showers due to decaying 
DM models are comparable to those for the annihilat- 
ing DM models. For > 1 TeV (m^ > 4 TeV), the 
shower rates for decaying DM models are larger than for 
the Kaluza-Klcin (Icptophilic where xX ~^ t^~^ f^~) model. 

The time required for the shower signal to be at 2a 
level for different values of the cone half angle centered 
at the Galactic center is presented in Fig. [15] for the 
Kaluza-Klein DM particle (solid line) and for Icptophilic 
model (dashed line) as example for annihilating and de- 
caying DM models. We find that the optimum time is for 
the cone half angle of 4° - 10°. For Icptophilic model, one 
could reach 2a signal within few months of observation 
for Omax = 10°. In case of Kaluza-Klein DM particle, 
one needs about 4 years of data taking for ^^,nax ~ 4°. 
For comparison, in Fig. [16] we show time needed for de- 
tection of 2(7 signal for the leptophilic model when DM 
particle decays (dashed line), gravitino two-body decay 
(solid line) and gravitino three-body decay (dotted line) 
models. The optimum angle for decaying DM models is 
about 50°. This is similar to the previously discussed 
muon signal presented in Fig. 13, but shower signals 
seem to be better for possible detection of the gravitino. 
This is because of better statistics with the increase in 
the cone half angle and also the inclusion of the con- 
tributions from the chargcd-currcnt and neutral-current 
interactions of the tau neutrinos and electron neutrinos 
with the detector medium. We note, however, that the 
20 years required for 9 > 50° for these parameters is not 
feasible in practice. 

In Fig. [TS] we present results for DM annihilation 
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FIG. 16: The time t versus cone half angle Smax centered 
at the Galactic center to reach 2a detection significance for 
the hadronic showers produced by the neutrinos from x 
f^^ l-i^ (dashed line) and from the two-body (sohd line) and 
three-body (dotted line) decay models of gravitino. We take 
El'^ = 50 GeV. 

cross section required for a given DM mass in order to 
reach 2a detection significance with shower events in five 
years of observation within Omax = 1°- From Fig. [T7] 
we note that the shower event rates decrease with 
for annihilating DM particles for fixed annihilation cross 
section, therefore the curves for DM annihilation cross 
section versus DM mass in Fig. [18] increase in order for 
reach 2a effect for a given m^. The parameter space 
above the curves defines the exclusion region and the 
leptophilic model seems to be more constrained than the 
Kaluza-Klein model. For the decaying DM models, our 
results for shower event rates obtained with fixed decay 
times (Fig. [T7| indicate that the 2a detection significance 
in five years curves for the decay time versus DM mass 
should increase with for < 1 TeV and becomes 
almost flat for niy, > 1 TeV. This is shown in Fig. [TO] In 
this case, the exclusion regions for the model parameters, 
decay time and the mass, at 2a level in five years are 
the regions below the curves. For the leptophilic models 
with parameters that satisfy Eq. (2), 2a signal detection 
would imply that the DM particle mass is 250 GeV in 
case of the annihilating DM, while for the case of the 
decaying leptophilic DM, the DM mass would be 3 TeV. 



V. CONCLUSIONS 

We have studied neutrino signals from DM annihilation 
and decay in the Galactic center assuming NFW profile 
as the DM density distribution in the Galaxy. We have 
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FIG. 17: Shower event rates as a function of DM mass for 
models considered in Table II. The boost factors for annihi- 
lating DM models and the decay times for the decaying DM 
models are fixed to the values in Table II. We take Smax ~ 50° 
and Elt = 50 GeV. 



FIG. 18: Annihilation cross section versus DM mass for 
shower events to reach 2a detection significance after five 
years of observation for the case of Kaluza-Klein (solid lines) 
and leptophilic (dashed lines) models. We take Smax ~ 50° 
and El'l = 50 GeV. 



considered models in which DM particle is a gravitino, 
a Kaluza-Klein particle and a particle in a leptophilic 
model. In case of the leptophilic model, we have consid- 
ered both the case of decaying and annihilating DM. For 
a gravitino, we have taken into account both two-body 
and three-body decay channels. For each DM model, we 
have calculated contained and upward muons and show- 
ers using the model parameters that were obtained by 
fitting the excesses in 7-ray and in the positron or elec- 
tron plus positron data from the observations of HESS, 
PAMELA and FERMI/LAT. We have used a range of 
cone half angles for the muon and shower events, and we 
have studied the dependence on the choice of the cone 
size and DM mass. 

Our results are summarized in Table I VII The specific 
models which are designed to account for the lepton ex- 
cesses, listed in Table II, are indicated by the italic entries 
(highlighted in red in the online version). In addition, in 
Table I VII we show the event rates and time required for 
a 2(7 observation for a range of masses and several val- 
ues of 6'max- Note that the event rates for different 0,„ax 
can be obtained by using rescaling of J-factors, for ex- 
ample Ji(50°)/Ji(10°) = 9 and Ji(10°)/Ji(l°) = 50 for 
decaying dark matter models; J2(50°)/ J2(10°) — 2.7 and 
J2(10°)/J2(l°) = 7.2 for annihilating dark matter mod- 
els. In Table IVlIl one can find the values of J(6'max)- 

For the leptophilic model {xX ^ M^y" ); for example 
with m^ = l TeV and B = 400, the muon flux due to DM 
annihilation dominates over muons produced by the at- 
mospheric neutrinos, for the muon energies in the range, 
200 GcV <Ef,< 950 GeV (50 GeV< <750 GeV) for 



the contained (upward) muon events with cone half angle 
^'max = 1°. The shower event rates for this same model 
never exceed the atmospheric neutrino induced shower 
rates, for 6'niax = 50°, in part due to the cone size. 

For = 800 GeV, B = 200 and 6'max = 1°, the 
contained muon flux due to annihilating Kaluza-Klein 
particles becomes larger than the background for i?^ > 
300 GcV and up to the kinematic cut-off i?^ = to^ = 800 
GcV. The muon flux is comparable to the background for 
muon energies 200 GeV < Ef^ < 400 GeV for the upward 
muons. 

With the model parameters that wc considered, the de- 
caying DM models do not produce the muon signal (up- 
ward or contained) that is above the atmospheric back- 
ground. The muon flux from the decaying DM would 
be comparable to those for the annihilating DM models 
only if the decay times were of the order ~ 10^^ sec and 
^ 10^^ sec for contained muon and upward muon events, 
respectively. In contrast to the muon case, we find that 
for a wide range of shower energies, the shower flux for 
decaying leptophilic particle (x — > A'^M ) is larger than 
for annihilating Kaluza-Klein particle. 

We have also calculated the total muon and shower 
event rates by folding the corresponding fluxes with the 
energy independent IceCube/DeepCore effective volume, 
i.e Vcs = 0.04(0.02) km'^ for the contained muon (shower) 
events. Due to its location IceCube detector is unable 
to study the upward muons produced by the neutrinos 
coming from the Galactic center. However, we have cal- 
culated the upward muon rates for the IceCube-type de- 
tector in the Northern hemisphere by folding the muon 
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FIG. 19: Decay time versus DM mass for both shower events 
to reach 2a detection significance after five years of obser- 
vation for the case of decaying DM models: gravitino three- 
body (solid lines), gravitino three-body (dashed lines) and lep- 
tophilic (dotted lines). We take 6'max = 50° and El'^ = 50 
GeV. 



fluxes with a muon effective area, which is assumed to be 

Aoff = 1 

Even if there is a significant signal to background ra- 
tio, low statistics may yield difficulties in confirming the 
presence of a DM signal via neutrinos. Thus, we have 
evaluated how many years it would take to observe 2a 
effect. Using our results for the muon and shower event 
rates, we have also obtained the contour exclusion plots 
in which we show the regions for the model parameter 
space for each DM model in the case of no signal detec- 
tion at 2(7 detection significance in five years. 

We find that the leptophihc model {xx A'^y" ) 
stronger constraints on DM annihilation cross section (or 
the boost factor) and than the case of the Kaluza- 
Klein particle. In terms of the constraint on the annihi- 
lation cross section, for the leptophihc DM model where 
the boost factor and DM mass are related by Eq. (2), 
after five years, the range of niy^ > 250 GeV would be 
excluded by upward muon events for Acs = 1 km^ and 
events within a cone half angle of ^max = 1°. A simi- 
lar limit is obtained from the shower rate as well. More 
generally, we find that the upward muon and the shower 
events arc more constraining than the contained muon 
ones. If there is no upward muon signal detected in 
five years, with ^max = 10°, for the decaying leptophihc 
model {x A'^/^ ) which satisfies the constraint given 
by Eq. (2), is constrained to be smaller than 3 TeV. 

In our calculations, we have taken the detector muon 
and shower energy thresholds to be 50 GeV. Changing 
the detector threshold energy to about 100 GeV does not 



affect our results significantly. However, decreasing the 
threshold energy down to ~ 10 GeV results in a larger at- 
mospheric background relative to the DM signal. Conse- 
quently, detecting a DM signal via muon or shower events 
with low detector thresholds 10 GeV) becomes more 
difficult. 

Increasing the cone half angle, 6'max about the Galac- 
tic center increases the DM signal, however it docs not 
necessarily improve the detection significance since the 
background signal due to atmospheric neutrinos is also 
enhanced. We have found the optimum cone half angles 
for all types of events in order to reach 2a detection level. 
For the annihilating DM models, we have found the op- 
timum angle to be a few degree for the muon events and 
about 10° for the shower events. In both cases, we have 
shown that there is a good chance of detecting both lep- 
tophihc and Kaluza-Klein particles in less than ten years 
for some DM masses. 

In the case of the decaying DM models, the optimum 
angle is about 50° for both muons and showers. For 
gravitino DM, signals could be detected in ten years only 
with shower events, while the decaying leptophihc parti- 
cle can be detected with upward muon events as well in 
few years. 

Our results with fixed annihilation cross section or de- 
cay time and variable DM mass can be used to predict 
neutrino signals for any other values for the parameters 
of the models we consider. For example, in the litera- 
ture, the HESS data is also explained by a hypothetical 
Kaluza-Klcin particle with mass = 10 TeV [3]. The 
boost factor that is required, in this case, is B = 1000. 
By rescaling the results in the previous sections, we find 
that the observation time to reach 2a detection signif- 
icance becomes 0.1 years for the upward muon events. 
For contained muons and showers, it is not feasable to 
detect 2a effect within reasonable time {t ^ 20 years). 
Increasing the boost factor by a factor of 5 and the DM 
mass by a factor of about 10 relative to the parameters of 
Table |TT] significantly improves the chance for detecting 
the Kaluza-Klein particle via upward muon events. 

Similar to the decaying leptophihc model, the grav- 
itino model ("03/2 l^l^v) with a gravitino mass of 3.3 
TeV and decay time r = 5 x 10^^ sec can also account 
for the FERMI data [13] ■ For a decaying DM particle, 
increasing the DM mass (for a fixed decay time) increases 
the event rates. Since the neutrino flux scales as ^ 
(see Eq. (|A2p ). decreasing the decay time also enhances 
the neutrino signals for a decaying DM particle. There- 
fore, combination of higher mass and shorter decay time 
should increase all the event rates that we have calculated 
for a lighter {rri^^^^ = 400 GeV) gravitino particle which 
has a longer decay time (r = 2.3 x 10^^ sec). For the 
model parameters, = 3.3 TeV and r = 5 x 10^^ sec 
we find that the 2a detection significance can be reached 
in 2.6 years via upward muons and in less than a year 
via the hadronic showers. In addition, for the contained 
muons the observation time decreases by two orders of 
magnitude relative to the value given in Table IIVI 
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9.2 


9.2 


9.2 


Br = 2.3 


iv;^*'(5o°) 


9.9 


50.9 


95.6 


139 


181 


364 


638 


844 


1010 


1150 




iV3h(50°) 


3.6 


7.66 


9.6 


10.74 


11.5 


13.17 


14.12 


14.46 


14.64 


14.74 




t7(10°) 


1 X 10* 


378 


107 


51 


30 


7.5 


2.5 


1.4 


1 


0.8 






2693 


101 


29 


14 


8 


2 


0.7 


0.4 


0.3 


0.2 




tsh{50°) 


210 


47 


30 


24 


21 


16 


14 


13 


13 


13 




iV^'(50°) 


2.13 


6.45 


8.43 


9.5 


10.2 


11.5 


12.2 


12.4 


12.5 


12.6 


= 2.9 


Af;^*'(50°) 


3.14 


29 


62.3 


97 


131 


286 


533 


728 


886 


1022 




iV.h(50°) 


1.95 


8.22 


12.09 


14.55 


16.2 


20.2 


22.45 


23.27 


23.68 


23.94 




^7(10°) 


1 X 10^ 


1 X 10^ 


252 


104 


57 


12 


3.5 


1.9 


1.3 


0.97 






2.6 X 10 


^ 316 


68 


28 


15 


3.2 


0.93 


0.5 


0.34 


0.26 




tsh{50°) 


709 


40 


19 


13 


11 


6.9 


5.5 


5.2 


5 


4.8 


-> ... 


iV^*(10°) 


14.2 


9.8 


7.2 


5.6 


4.6 


2.4 


1.25 


0.84 


0.63 


0.51 


B = 200 


Af^f(10°) 


86.1 


131 


140 


130 


128 


124 


108 


92 


81 


72 




A^.ft(10°) 


11 


9 


7 


5.7 


4.8 


2.6 


1.4 


0.9 


0.7 


0.6 




^7(1°) 


1.27 


0.63 


0.54 


0.65 


0.66 


0.7 


0.87 


1.14 


1.42 


1.72 




C(io°) 


1.55 


0.68 


0.57 


0.71 


0.72 


0.76 


1.0 


1.36 


1.76 


2.2 




C(50°) 


5.1 


2.2 


1.84 


2.29 


2.3 


2.44 


3.2 


4.5 


5.8 


7.2 






3.4 


4.4 


5.9 


7.7 


9.6 


22 


61 


116 


189 


280 




t.;^(10°) 


1.3 


1.9 


2.9 


4.3 


5.8 


18 


64 


136 


237 


364 




tsh{50°) 


3.3 


5 


8 


12 


16.3 


57 


204 


445 


777 


1202 


XX -> m'^m" 




40.19 


29.58 


22.01 


17.39 


14.3 


7.59 


3.90 


2.63 


1.98 


1.59 


B = 400 


n;:p{io°) 


144 


241 


273 


283 


320 


266 


221 


190 


167 


151 






51.4 


45.6 


36.4 


30 


25 


14 


7.4 


5 


3.8 


3 






1.11 


1.68 


2.55 


3.61 


4 


13.64 


44 


92 


156 


238 




i?(10°) 


0.66 


1.18 


2.06 


3.24 


4.7 


16.31 


61 


133 


234 


364 




t?(50°) 


1.93 


3.55 


6.38 


10.2 


15 


53 


201 


444 


781 


1213 






0.54 


0.24 


0.2 


0.18 


0.14 


0.21 


0.28 


0.35 


0.43 


0.50 




t7(10°) 


0.47 


0.21 


0.16 


0.15 


0.12 


0.17 


0.25 


0.33 


0.42 


0.52 




C(50°) 


1.83 


0.65 


0.51 


0.47 


0.37 


0.54 


0.78 


1.1 


1.35 


1.7 






0.63 


0.72 


0.91 


1.12 


1.37 


2.58 


5.5 


9 


13 


18 




f«h(10°) 


0.12 


0.14 


0.2 


0.26 


0.34 


0.87 


2.63 


5.34 


9 


13.6 




t.h(50°) 


0.18 


0.22 


0.33 


0.48 


0.7 


2.1 


7.2 


15.5 


27 


42 


Atmospheric 




2.28(1°) 






227.5(10°) 




5347(50°) 








28(1°) 






2794(10°) 




65668(50°) 








0.3(1°) 






28.8(10°) 




676(50° 


) 



TABLE VI: Summary of the resuhs for the event rates and the time that it takes to reach 2a effect for different values of 
and 6'max, where the threshold enery is taken as 50 GeV for both muon and shower events. In this table we do not include 
results which have f > 15 years for all m^. Results for the specific choice of the parameters in each model corresponding to 
fitting PAMELA, FERMI/LAT and HESS data are presented in red italic fonts. 
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The dependence of the signal on the cone half angle is 
different for the annihilating DM particles than for the 
decaying DM particles. We have demonstrated this by 
choosing wedges between ^max — 1° and 0niax centered 
at the Galactic center and calculating (J„)nAri, which 
defines an overall normalization for the neutrino signals, 
for different Omax for both annihilating and decaying DM 
particles. The angular wedges can be used to rescale 
event rates as well. Our results indicate that the J fac- 
tor for annihilating DM particle decreases sharply with 
^max whereas for the case of decaying DM particle, for a 
wide range of ^max, the J factor has a weak ^max depen- 
dence. In determining the nature of the DM (annihilating 
or decaying), the directional dependence of the neutrino 
signals gives valuable information. 
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Appendix A: Neutrino Flux 

The neutrino flux observed at the Earth due to the DM 
annihilation or decay in the galaxy for a given neutrino 
flavor can be written as 1231 



[av) 



^ \ F 



E^f^^ I <•^2)^2A^^ (Al) 



for the case of annihilating DM, and 



dE^ 



RoPo 



dNl 
dE,, 



{Ji)n^n (A2) 



for the case of decaying DM where dN^ /dE^ is the neu- 
trino spectrum for a given annihilation or decay channel 
F with branching fraction Bp, B is the boost factor, r 
is the decay time, Ro is the distance of the solar sys- 
tem from the Galactic center and po is the local den- 
sity near the solar system. The neutrino energy spec- 
trum, dN^/dEij for different channels can be found in 
Appendix B. In the above equations, the dimensionless 
quantity (J„)n is defined as p^.[2^[26j 



where p{l) is the DM density, l{9) is the distance between 
the source and the Earth in the direction of 9 which is the 
cone half angle from the Galactic center and the integral 
is over the line of sight (l.o.s) within a solid angle Ail = 
27r(l — cosdmax), Centered in the Galactic center. 

In our calculations, we take the DM annihilation cross 
section to have the typical thermal relic value {av) = 
3 X 10~^^ cm^s~^ and we use the Navarro- Frenk- White 
(NEW) DM density profile [H 



p(r) = 



{r/Rs){l + r/R,Y 



(A4) 



where ps and Rs are the parameters which vary from halo 
to halo. In our calculations we set ps = 0.2589 GeV/cm'^ 
and Rs = 20 kpc so that the DM density in the vicinity 
of the solar system (r = Ro = 8.5 kpc) takes the typical 
value p{Ro) = 0.3 GcV/cm'^ ^^l- Using these definitions, 
we can write r = yj R^-\- P — 2Rol cos and the upper 
limit for the I integral in Eq. ()A3|) can be obtained as 

^max — Ro COS 6 + \J Rl — R^ SVC? 9 for a given 9. 

For the NFW profile, some values for (J2)oAr2 and 
(Ji)oAf2 are summarized in Table IVnl 





0.1° 


1° 


5° 


25° 


50° 


70° 


90° 




0.14 


1.35 


5.94 


19.68 


27.75 


31.73 


33.42 


(Ji)nAJ7 


0.00027 


0.018 


0.30 


3.69 


8.79 


12.24 


14.90 



TABLE VII: The values of J factors for NFW profile for 
Bma^ = 0. 1° , 1° , 5° , 25° , 50° , 70° , 90° . 

In Fig. [20] we show (J„)f2Ari factors as a function of 
6'max for DM annihilation {n = 2) and DM decay {n ~ \) 
evaluated for a cone wedge between 0max — 1° and ^max 
around the Galactic center. 



Appendix B: Neutrino Spectra 

In this study, we have studied the model dependent 
neutrino signals from the annihilation/decay of the DM 
particles that reside in the galaxy. We assumed that 
the DM particles are non-relativistic so that their total 
energy is close to their rest energy {E^ ~ "^x)- 
neutrinos with energy E^, can be produced from the an- 
nihilation/decay of the DM or from the decay of quarks, 
charged leptons and gauge bosons which are produced by 
the annihilation/decay of the DM. In our calculations, we 
have used the standard unpolarized decay distributions 
which, in general, take one of the following forms 



dN_ 
dx 



25/ (3x2 _ 2x3) 



(Bl) 



{Jn)n — 



r dn 


r di{9) 




1 All J 


l.o.s Ro 


if) 



(A3) 



dN 
dx 



12%(a;^ -a;^) 



(B2) 
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e (degrees) 

max ^ ^ 



and 

T.z = —Tz. (B7) 
We take Tz = 2.5GeV. 

2. Neutrino spectrum from r^, fi^ decay channels 

In this section, we present the spectrum from ^ and 
T decays. The spectra for other neutrino flavors can be 
deduced from these results. For example, the Ve spectrum 
from ^ decay is identical to the i'^ spectrum from r decay 
and the Vr spectrum from r decay to that of from fi 
decay. For these three-body decays, the j/^ spectrum 
from T decay is given by Eq. (jB2[) and from fi decay by 
Eq. (|B1[) in the fr'amc where the decaying particle is at 
rest. Then, after boosting these results by using Eq. (|B4p 
for charged leptons produced via the DM annihilation or 
decay, we obtain 



FIG. 20: J factor values both for annihilating (dashed line) 
and decaying (solid line) DM models for a wedge between 
5max — 1° and Smax about the Galactic center as a function of 



dN 
dx 



BfS{x - 1) 



(B3) 



in the rest frame of the decaying particle where x = 
2E^/md, rrid is the mass of the decaying particle and Bf 
is the decay branching fraction for a given decay chan- 
nel. Once the distribution in the rest frame is known, the 
neutrino energy spectrum from a decaying particle with 
velocity Pd and energy Ed = ^drrid is given by [iij 



dN, 
dE, 



de 1 UN 

e Idfid V de 



(B4) 



where E^^E,-f^\l±l3d)-^. 



dN^ 



El ( 3 "^^^ 



(B8) 



for the case of anni- 



where x ~ < 1 and Ei ~ 
hilation or Ei = for the case of decay. The decay 

branching fraction, Bf = 0.18(1) for r(/i) decay. In addi- 
tion to the three-body decays, r can also decay into Vr via 
T — > Vt-M or t — >■ i^T-X where M = 7r,p, ai mesons and 
X indicates hadrons. The neutrino spectra from these 
channels are given as [iHl 



dE„. 



B 



J-— when ^ < rM for mesons 

l-M Er (B9) 

when 1^ < 0.3 for hadrons , 



0.3St 



where rjv/ = 1— m^^/m^ with tom and mj- being the mass 
of the meson M and the r Icpton, respectively. Here, 
Bf = 0.12,0.26,0.13 and tm = 0.99,0.81,0.52 for M = 
TT, p, oi respectively. 



1. Neutrino spectrum from x ^ Zi' decay channel 



Neutrino spectrum from 6(6) and c(c) decay 
channels 



In our calculations, for the chann el y — >■ Zv we use the 
Brcit-Wigncr distribution given as [23| 



1 



dN, ^ 

dE, - {El-El,f + El,Tl, 

dE 



where the distribution peaks at 



E„ 



, (B5) 



(B6) 



The b and c quarks hadronize before they decay into 
neutrinos. The hadronization effect is taken into account 
by scaling the initial quark energy, Ein , in the form Ef ~ 
ZfEjn, where f = b, c, Zf = 0.73(0.58) for b (c) quarks 
[451 1 and Em = for an annihilating DM particle or 
Ein = for a decaying DM particle with mass m^. 

The neutrino spectrum from the decay of /= 6, 6, c or 
c from XX ~^ ff can also be ap pro ximated by the second 
equation in Eq. (|B8|) (see also [H, ) , i.e. 



dN, 
dE, 



2B 



(1 - 3x^ + 2x^) for a; < 1 



Ef 

otherwise 



(BIO) 
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where x = and 



(Ef , Bf) 



(O.lSEin , 0.103) b channel 
{0.58Ein , 0.13) c channel. 



(Bll) 



4. and Z decay channels 

In the and Z decay channels, the neutrino spec- 
trum from the decaying particle with velocity, /3b, and 
energy, Eb, can be obtained by using Eqs. (|B3p and 
(El, i.e, 



^ for ^il-h)<E.<^{l + M 



dE,, 



otherwise , 



where 



(1 , 0.105) W channel, 
(2 , 0.067) Z channel. 



(B12) 



(B13) 



5. Neutrino spectrum from t{t) decay channel 

The top quark decays into a W boson and a b quark 
{t — > Wb) with a branching fraction close to unity. Thus, 
the sum of neutrino spectra of W and b channels gives 
the required spectrum, i.e, 



where Bw = 0.105,5^ 

E^ ~ ^ -u « ^-l 



0.103,eT = EwilTPw)/2, 
E^j-\l ± I3t)-^ and E^ = 0.73Eb with Ew, Pw 
and Eb being equal to their values in the top-quark rest 
frame, i.e. 



Eb 
Ew 



2mt 
rnf + m 



W 



2mt 



w 



Eb 
Ew 



(B18) 



The function Db is given by 
1 



Dblm, n] = - [gfm^ - n^) - 4(m^ - n^) + 6ln ( —) 

(B19) 



6. Neutrino spectrum from ipz/2 — > I I v decay 
channel 



In the zero mass limit, the primary lepton (l^ or or 
v) spectrum from the decay channel, '03/2 l^l^iy, can 
be approximated to be 



dN, 



60 



dE, 



l{u) 



-x^{l — x) where x 



3/2 



2E 



< 1 , 



(B20) 

by using the results in [48|. In order to obtain the spec- 
trum for the secondary neutrinos produced from the pri- 
mary charged lepton decays, one can use 



dN^ 
dE,, 



[dN, 



w+w- 



fdN^ 

\ d'B„ / 



(B14) 



Then, boosting this expression yields the neutrino spec- 
trum from top quarks moving with velocity /3f and energy 

Et = 7t"it m, 



dN^ _ UN,,'- 
dEn \ dE,, y 



[dE, 



(B15) 



/ b 



where 



dN„ 
dE,, 



B 



w 



-h 



min(£'+, e+) 



y / w 



2^tPtEwPw max(i?_,e_) 
if 7t(l - /3t)e_ <E,< 7t(l + /3t)e+ 
otherwise , (B16) 



and 



dN,, 



dE, 



Bb 



/ b 



2-ftPtEd 
if E, < 7t(l 
otherwise 



Db[E-/Ed,mm{l,E+/Ed)] 
Pt)Ed 



(B17) 



dN„ 
dE„ 



dEi 



1 dNi 
Ni dEi 



dN, 
dE„ 



(B21) 



Here, is the primary charged lepton spectrum given 
by Eq. (|B20|) and the spectra, I ) are given by 

Eqs. (jB8|) and (jB9p . Then, the secondary i/^ spectrum 
is derived to be 



fdN„ 



\dE„^ 



5B 



f 



from the primary /i decays and 



(1 - 6x2 + g^3 _ 3^4^ ^ (-g22) 



dE,^, 



6B 



f 



(l-10a;2 + 20a;^-15a;* -1-4x5) ^ (323) 

"^V'3/2 

from the primary r decays, requiring that x < 1 in each 
case where x = 2E,j/m^,j^^^. Finally, for the Vr spec- 
trum from the primary r decays accompanied with the 
meson/hadron production, we find 



dN„^ 



3Bf 



1 - 



with the requirement tm > x. 



M 



(B24) 
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Appendix C: Parametrizations 

1. Contained and Upward Muons 

Our results for the contained muon fluxes given in Fig. 
4 can be parametrized as 



B_\ / (J2)A0 \ gx) 
lOOy V 1.35 J {m^/TeVy 



Vl026 



sec 



0.018 



(mx/TeV) 



for annihilation and decay processes, respectively. The 
upward muon fluxes presented in Fig. 6 can be 
parametrized as 



^7 



dE^ 
dE,, 



100 



1 



c{x)ax) 



V 1026s 



1.35 J (m^/TeV) (i + 1.5^ 
{Ji)An\ C{x)gx) 



0.018 / 



(1 



1 c,-^V 



(C2) 



where x = Ep^/m^ for DM annihilation and x — 2ii'^/mj^ 
for the DM decay. The functions ^(x) and C{x) are fitted 
to our results for muon fluxes and parametrizes as 



ai + 026 

„4 



+a6a;^ + a-j ln(x) 
C{x) = 1 + ag,x , 



+ + a^x^ + 052;'^ + 



(C3) 



with the best fit values for the parameters given in Tables 
(jVIIip and (jlXp for the contained and upward muons, 
respectively. 





V'3/2 


^"3/2 


X 




XX 




(three-body) 


(two-body) 




-¥ ■ ■ ■ 




ai 


16.651 


-32.032 


0.0554 


2.5637 


0.029 


a2 


-16.642 


32.04 


-0.0472 


-2.5470 


-0.0012 


a-i 


-16.640 


32.03 


-0.0472 


-2.5546 


-0.0012 


a4, 


8.242 


-15.93 


-0.00934 


1.2381 


-0.11 


0,5 


-2.553 


5.022 


0.0254 


-0.3780 


0.11 


ae 


0.423 


-0.87 


-0.0069 


0.06929 


-0.027 


a7(10-'') 








-0.0765 


215.8 


-0.04 



TABLE VIII: The best fit parameter values for the contained 
muon fluxes for different DM models. 

The background muon flux from atmospheric neutrinos 
can be written in a parametric form as 



dE„ 



= 3186.2 



M / ATM.avE 



vGcv; 



An 

10-3sr 



(C4) 

in units of GeV ^km ^yr ^ for the contained muon 
events and 





^3/2 


^3/2 


X 




XX 




(three-body) 


(two-body) 


fi+fi- 


■ ■ ■ 




ai 


1.251 


1.292 


-0.7716 


-3.4381 


-4.60 


a2 


-1.239 


-1.285 


0.781 


3.4847 


4.67 


as 


-1.266 


-1.296 


0.747 


3.3390 


4.44 


04 


0.638 


0.641 


-0.341 


-1.5573 


-2.0 


as 


-0.203 


-0.196 


0.0918 


0.4439 


0.526 


aa 


0.0364 


0.032 


-0.0132 


-0.0694 


-0.073 


^7(10-'^) 








0.45 


41.47 


3.6 


as 


0.3 


0.3 


1.5 


1.2 


1.5 



TABLE IX: Same as Table (|VIII|) but for the upward muon 
fluxes. 



up 



= 89.0 



(. 



E„ 



-1.475 



An 



(C5) 



in units of GeV~^km"^yr^^ for the upward muon events. 



2. Showers 

Our results for shower flux presented in Fig. 14 can be 
parametrized as 



dEsh 
d4>sh 



B 
100 
r 



dE. 



sh 



(l026sec) ( 



MAn 

27.75 ) (mx/100GeV)2 
-i/(Ji)Ar!\ f(x) 



8.79 J (mx/GeV) 



(C6) 



for DM annihilation and DM decay processes, respec- 
tively. 





V'3/2 


V'3/2 


X 


5(1)5(1) 


XX 




(three-body) 


(two-body) 




. . . 




ai 


-6.335 


-14.129 


-22.867 


-8.765 


-32.43 


a2 


6.355 


14.139 


22.898 


8.78 


32.48 


as 


6.274 


14.100 


22.765 


8.72 


32.29 


a4 


-2.991 


-6.952 


-11.134 


-4.25 


-15.79 


as 


0.813 


2.127 


3.325 


1.25 


4.71 


ae 


-0.099 


-0.347 


-0.5125 


-0.189 


-0.73 


a7(10-=*) 


-0.472 















TABLE X: The best flt parameter values for the shower fluxes 
for different DM models. 

The background shower flux from atmospheric neutri- 
nos can be written in a parametric form as 



[>sh 



3.21 X 10' 



-2.155 



ATM,avg \GcV/ 

in units of GeV^^km~'^yr~^. 



An \ 

2.24sr ) 
(C7) 
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